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Abstract 

Range size variation in closely related species suggests different responses to 
biotic and abiotic heterogeneity across large geographic regions. Species turn- 
over generates a wide spectrum of species assemblages, resulting in different 
competition intensities among taxa, creating restrictions as important as envi- 
ronmental constraints. We chose to adopt the widely used phylogenetic related- 
ness (NRI) measurement to define a metric that depicts competition strength 
(via phylogenetic similarity), which one focal species confronts in its environ- 
ment. This new approach (NRIfo^ai) measures the potential of the community 
structure effect over performance of a single species. We chose two ecologically 
similar Peucaea sparrows, which co-occur and have highly dissimilar range size 
to test whether the population response to competition intensity is different 
between species. We analyzed the correlation between both Peucaea species 
population sizes and NRIfoj-ai using data from point counts. Results indicated 
that the widespread species population size was not associated with NRIfocah 
whereas the population of restricted-sized species exhibited a negative relation- 
ship with competition intensity. Consequently, a species' sensitivity to competi- 
tion might be a limiting factor to range expansion, which provides new insights 
into geographic range analysis and community ecology. 



doi: 10.1002/ece3.959 



Introduction 

Biologic interactions are fundamental to the understand- 
ing of ecological patterns and processes. Webb et al. 
(2002) set the framework of phylogenetic community 
ecology to elucidate the effect of ecological interactions 
depicted via phylogenetic structure, that is, the degree of 
relatedness among species forming an assemblage. This 
approach is strongly supported by the tested observation 
of the inheritance of niche-related traits from ancestors, 
termed niche conservatism (Wiens et al. 2010; for birds 
see Lovette and Hochachka 2006), which reinforces the 
concept of a link between phylogenetic distance and eco- 
logical similarity reflected in behavioral (Houle 1997) and 
life-history traits (Burns and Strauss 2011). Consequently, 
demonstrating statistical support for phylogenetic and 



community assemblage patterns (i.e., clustering or even- 
ness) suggests a potential process exists structuring species 
assemblages (Webb et al. 2002). Questions addressed 
under this method have resolved several issues related to 
ecology (ecosystem stability Cadotte et al. 2012; species 
lost and climate change, Willis et al. 2008), providing 
useful insights into the internal structure of phylogenetic 
and ecological relationships (Lovette and Hochachka 
2006, Gomez et al. 2010). For example, Graham et al. 
(2009) demonstrated that biologic interactions among 
hummingbird species represented the leading factor in 
tropical lowland community assembly processes, even 
more important than environmental variation (i.e., habi- 
tat filtering). However, a bond must exist between com- 
munity phylogenetic structure and the performance of 
each species, a fact that is largely overlooked. The fitness 
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of particular species might consequently be explained by 
the composition of the community in which the species is 
found. Furthermore, Ricklefs (2004) reported that com- 
munity level processes generated population level changes, 
driving current ecological patterns. Yang et al. (2013) 
evaluated the phylogenetic diversity (PD) of assemblages 
surrounding target individuals, showing that most indi- 
viduals had a neutral tendency regarding the PD of adja- 
cent plots. However, by exploring the phylogenetic 
distance of each individual in a focal species, the direct 
effect of the surrounding community structure may be 
measured. For instance, Jiang et al. (2010) designed an 
experiment to assess the success of an invader species in 
bacterial communities, based on phylogenetic distance 
between invader and native species. The authors found a 
positive relationship between phylogenetic distance and 
the probability of a species becoming established. 

The central role of biotic interactions is considered 
contingent on a species geographic range. For example. 
Brown et al. (1996) reported that biotic interactions tend 
to limit the distribution and abundance of species at 
lower latitudes. Differences in a species geographic range 
size do not just suggest variation in response to environ- 
mental variables or niche breadth (Gaston and Spicer 
2001), but can reflect a species response to biologic inter- 
actions, which were illustrated in classic experiments by 
ConneU (1983), and more recent studies by Bullock et al. 
(2000). However, previous studies reported the influence 
of biologic interactions in two-species systems distribu- 
tions, without evaluating the effects of the entire commu- 
nity. Recently, Villalobos et al. (2013) introduced a novel 
and interesting approach, in which the phylogenetic struc- 
ture of species co-occurrence of a focal species is used to 
study broad coexistence patterns. 

We hypothesize that the sensitivity of species confront- 
ing negative interactions is reflected in the species popula- 
tion attributes: species inhabiting different assemblages, 
consequently experience different levels of competition 
throughout the species geographic range, which results in 
different rates of change in a species overall fitness. Gas- 
ton (2009) indicated that population size is the outcome 
of several population structure attributes (e.g., levels of 
births, deaths, and migration). For instance, population 
density has been linked to species richness, producing 
higher density in areas where richness is low (i.e., density 
compensation; MacArthur et al. 1972), which is a pattern 
that was first described for island systems compared with 
mainland systems. Under these conditions, it is expected 
that populations of a species that occupy large geographic 
ranges (widespread) are not as influenced by co-occurring 
species with which they compete; alternatively, popula- 
tions of species exhibiting restricted geographic ranges are 
more affected by increased potential competition. Support 
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for range size heritability (Waldron 2007; but also see 
Webb and Gaston 2005) facilitates the expectation that 
related species would have similar range size; however, 
differences in the range of closely related species might 
serve as a viable system to test whether this dissimilarity 
is provided by a differential response to competition. The 
present study included two components: (1) we modified 
a widely used metric of phylogenetic structure (NRI; 
Webb et al. 2002) to center the attention toward a focal 
species and (2) two sympatric species in the genus 
Peucaea (Emberizidae), which exhibit very dissimilar 
range sizes, were used to evaluate whether a population 
size response to potential competition (through the modi- 
fied metric) differs between the two sparrow species. 

Methods 

Peucaea sparrows and f ieldwork 

We conducted this study in southeast Mexico, in a region 
called the Tehuantepec Isthmus (Huidobro et al. 2006). 
This region is located in the narrowest stretch of land 
between the Gulf of Mexico and the Pacific Ocean, which 
is represented by the municipality of Juchitan, state of 
Oaxaca. Two sparrow species co-occur in the area: 
Peucaea ruficauda and P. sumichrasti. The former is a 
more widespread species (2.6 x 10^ km^; see next sec- 
tion) compared with P. sumichrasti, which is one of the 
most range-restricted avian species in Mexico 
(9.7 X 10^^ km^; Wolf 1977). This species is endemic to 
the Tehuantepec Isthmus, and its entire range overlaps 
with P. ruficauda. Both sparrow species are common birds 
in the region and are similar in shape, behavior, and ecol- 
ogy (Wolf 1977). We selected 17 monitoring sites across 
the region, which had comparable vegetation type, human 
perturbation, and size, covering almost the entire geo- 
graphic range of P. sumichrasti. The species composition 
(land birds) of each assemblage and Peucaea sparrow 
abundance were described by sampling 24 fixed-radius 
point counts separated by 200 m to avoid double count- 
ing individuals at each monitoring site. Each point count 
was sampled eight times by the same team of observers 
for one year, with each monitoring site being visited every 
6 weeks on average. 

Phylogeny and geographic range size 

We conducted a phylogenetic reconstruction, which 
included all land bird species we observed during our 
fieldwork. We queried the GenBank database (NCBI, 
September-October 2011) for the mitochondrial COI 
gene (Hebert et al. 2003; Alif et al. 2011; Appendix Al), 
representing each species we identified from our 17 
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monitoring sites. Sequences of the species not deposited 
in GenBank were replaced with an available sequence 
from the closest relative, according to previous studies 
reported in the literature. Sequence alignments were made 
in Clustal W (Larkin et al. 2007), a BioEdit 7.0.9.0 (Hall 
1999) accessory tool. The reconstruction was performed 
by Bayesian method using BEAST vl.7.5 (Drummond 
and Rambaut 2007). Nucleotide substitution model 
employed was HKY+G+I being identified as the more 
adequate by JModelTest 0.1.1. The Yule speciation process 
was set to model the tree prior. In order to calibrate the 
root node of the tree, we used the date from Jetz et al. 
(2012) for the divergence of Anatidae (lognormal mean 
4.27, SD 1, zero offset 0). One independent 5,000,000 
generation run was performed sampling at every 1000 
generations. The outcome was analyzed in TreeAnnotator 
vl.7.5 discarding the 10% of trees and visualized in 
Mesquite v2.74 (Maddison and Maddison 2011) and is 
available through Figshare (http://dx.doi.org/10.6084/m9. 
figshare.865723). 

We looked for range size phylogenetic signal as an 
Emberizidae family trait by conducting a randomization 
test (999 randomizations; Blomberg et al. 2003) using the 
comprehensive phylogeny of Emberizidae reported by 
DaCosta et al. (2009), and the range size of 49 sparrow 
species. The species distributions area calculations were 
generated in ArcGIS 10 and shapefiles from the Nature- 
Serve compilation (Ridgely et al. 2003) using the appro- 
priate geographic coordinates regarding species 
distributions in North or Central America. Phylogenetic 
signal is detected when random distributions exhibit sig- 
nificant differences from observed values. Although simi- 
lar geographic range sizes have been observed among 
close relatives in birds (Waldron 2007), we decided to test 
this hypothesis in the Emberizidae family because the out- 
come was scale dependent. 

Data analysis 

The traditional metric of the phylogenetic community 
(Net Relatedness Index -NRI-) is a standardized measure 
of the mean pairwise phylogenetic distance (MPPD, 
Webb et al. 2002), which is the phylogenetic distance 
among all possible pairs of species within a community. 
We modified the MPPD, so that our metric (MPPDfocai) 
did not reflect the distance among all pairs, rather it mea- 
sured the distance from the focal species and each species 
included in the assemblage, and averaged those values 
(Fig. 1). The calculation was performed using the cophe- 
netic distance between focal species and all others mem- 
bers of the community, from which the mean was 
calculated. The species abundances are included in the 
weighed version of this metric by calculating the weighed 
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Figure 1. The NRIfocai calculation using a phylogenetic tree and its 
phylogenetic distance matrix. Notice that the difference with 
traditional NRI appears in the MPPD calculation, where the traditional 
is estimated by averaging the distances among all the possible pairs, 
while MPPDfocai is estimated by averaging the distances among the 
focal species (species C in this example) and each species included in 
the assemblage. 

mean instead of the arithmetic mean; in this way, the 
metric reflects the real composition of the community. 
Like the NRI, we constructed a standardization named 
NRIfocai. in which the observed MPPDfocai values were 
compared with null distributions that were generated by 
creating communities of identical size by random draws 
from species pool (Kraft et al. 2007). In a single value, 
the NRIfocai describes focal species relatedness and the set 
of species that co-occur with the focal species. In this 
way, NRIfocai niay be defined as the phylogenetic (ecologi- 
cal) similarity of Peucaea with all other observed bird 
species and serves as a measure of potential competition 
(Fig. 1). 

We calculated NRIfocai (both weighed and non- 
weighed) for P. ruficauda and P. sumichrasti for each 
assemblage and performed a Pearson's product-moment 
correlation between NRIfocai values and sparrow abun- 
dance. Peucaea abundance was normalized by root square 
transformation (Sokal and Rohlf 1987). The modified 
metric was based on the "mpd" function implemented in 
the PICANTE package (Kembel et al. 2010; in R 2.15.1 
R Development Core Team 2010). The relationship of 
NRIfocai and NRI was calculated through 100 simulated 
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communities, measuring both the phylogenetic metrics of 
each assemblage and calculating the correlation between 
them. 

Results and Discussion 

The geographic range among Emberizidae sparrows repre- 
sents a trait with significant phylogenetic signal 
(P < 0.001), indicating that similar species have a similar 
range size (Fig. 2). However, the difference in the geo- 
graphic range size of sympatric and closely related 
Peucaea sparrows is noticeably large (2.5 x 10^ km^). 
This range size difference is of interest because of the eco- 
logical similarity between the sparrow species and the 
apparent lack of any form of geographic barrier; hence, 
this scenario may be explained by biologic interactions. 
We found that the widespread species {P. ruficauda) was 
more abundant compared with P. sumichrasti in 80% of 
the species assemblages. The relationship between range 
size and local abundance has been well established in 
several taxa and was demonstrated in passerine birds 
(Bock and Ricklefs 1983). Several mechanisms have been 



proposed to determine the positive correlation between 
range size and local abundance (revised and discussed by 
Gaston et al. 1997); for example, Hoh et al. (1997) con- 
ferred special importance to among-species differences 
based on species differential responses to density-indepen- 
dent factors influencing population attributes. Here, 
results emphasized that biologic interactions were an 
important element in this process. For the resident species 
assemblages, the number of P. sumichrasti individuals 
showed a negative correlation with NRIfoj-ai (r = —0.592, 
P < 0.05; Fig. 3) when considering all species abundances 
(weighed NRIfocai)- This suggests that competition 
strength (implied by increasing phylogenetic similarity) 
limits abundance in P. sumichrasti. The competition 
strength approached by the phylogenetic similarity leads 
the population volume in this restricted-sized species. 
Consequently, the community structure (composition of 
the species assemblage) is an important factor affecting 
one species population size. This result is analogous with 
that reported by Jiang et al. (2010), who found a positive 
relationship between the phylogenetic distance of invaders 
to bacterial communities and invader abundance. While 



Figure 2. Geographic range size among 
Emberizidae sparrows sinowed a pliylogenetic 
signal suggesting tliat closely related species 
have similar range sizes. Color indicates the 
size of geographic range (in thousands of 
square kilometers) for each species. 
White = 1-100; light gray = 101-200; dark 
gray = 201-500; slate gray = 501-1000; black 
<1000. 
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Figure 3. Scatter plot for correlation analysis between NRIfocai and 
Peucaea abundances (root square transformation). Tine species 
abundances are included in tine weiglied version of NRIfocai by 
calculating tine weiglned mean instead of tine arithmetic mean. 

Table 1. Correlation test results between Peucaea sparrow popula- 
tion sizes and our "focal" version of Net Relatedness Index (NRIfocai) 
from 17 bird assemblages. The relative abundances from each entire 
assemblage were used to weigh the NRIfocai. Note the negative rela- 
tionships represent an increase between population size and phyloge- 
netic dissimilarity. 





Non-weighed 

NRIfocai 


Weighed 

NRIfocai 






r 


P-value 


r 


P-value 


Peucaea ruficauda 
Peucaea sumichrasti 


-0.309 
-0.004 


0.227 
0.987 


0.2 
-0.592 


0.439 
0.012 



the context of this preceding study differs to ours, it is 
worth noting the similar outcomes, even for very different 
organisms and differently sized communities. However, 
we did not find evidence of a relationship between wide- 
spread P. ruficauda abundance and NRIfocai (Table 1). 
Dissimilarity in range size among ecologically similar spe- 
cies might be related to differences in sensitivity to com- 
petition among species. This difference might be linked to 
the categories defined by Yang et al. (2013) in a lightly 
similar context, who grouped phylogenetic diversity into 
neutral, repeUer, and accumulators species. Sensitivity to 
competition as a trait has the potential to limit species 
distributions across large geographic regions, because a 
large area exhibits high species turnover; consequently, 
several communities may form levels of competition so 
great that a sensitive species could not compete success- 
fully. Consequently, regions at low latitudes with high 
species turnover can be defiant to susceptible species 
producing small-sized range species and contributing 



(besides other factors; Arita et al. 2005) to the observed 
and proved pattern of Rapoport's rule. Yang et al. (2013) 
identified only a few phylogenetic diversity repellers, 
which might indicate that few sensitive species produce a 
reduced number of species with small-sized ranges and 
elevated numbers of species with large-sized ranges. How- 
ever, the opposite pattern is obtained in reality; whereby, 
many species have small- to moderate-sized ranges, while 
only a few have very large ranges (Brown et al. 1996). 
Thus, the complex interaction between biologic interac- 
tions and environmental variation is reaffirmed. 

The internal distributional range structure (sensu 
Brown et al. 1996), and the variability in species abun- 
dance throughout a species distribution (Brown et al. 
1995), might provide a way of elucidating the response of 
species to potential competition across different assem- 
blages. In even a relatively small area (Tehuantepec Isth- 
mus), the observed sites showed enough species turnover 
(see Appendix A2) to depict dissimilar phylogenetically 
structured assemblages, which suggested a gradient of 
potential competition (Fig. 4). Moreover, the geographic 
boundaries of sensitive species might have been estab- 
lished by communities in which ecological (i.e., phyloge- 
netic) similarity exceeds the level of similarity (i.e., 
competition) that a focal species could withstand. In 
other words, certain assemblages function as "stakes," 
limiting species geographic ranges. These stakes acted as 
biologic barriers setting a threshold that could not be tres- 
passed by sensitive species. For instance, the P. sumichrasti 
population size reached low levels when the community 
phylogenetic similarity was high, preventing the dispersal 
of individuals to new areas. Although some individuals 
might cross these stakes, their numbers might be reduced, 
due to their being unable to establish stable populations. 
Areas where phylogenetic similarity is low would facilitate 
the free transit of sensitive species, until another type of 
barrier (i.e., environmental or biologic) is reached. Even 
though our monitoring sites were distributed throughout 
the geographic range of the restricted species, it would 
have been more accurate to evaluate our hypothesis by 
identifying the exact boundary of the range to assess the 
presence of overwhelming competition to which sensitive 
species are subjected. These boundaries are not distinct 
for birds; however, it might be viable to apply such stud- 
ies that exhibit lower levels of movement to improve the 
experimental design. 

In our study, the specific assemblage composition alone 
was not enough to explain the variability in population 
size among sites; the positive relationship between 
P. sumichrasti abundance and the nonweighted NRIfocai 
was not statistically supported (P > 0.05). Communities 
might appear similar to each other when species' 
abundances are neglected. In the Appendix (A2), we show 
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Figure 4. Reconstruction of a phylogeny of the species included in the study and their presence at the 17 monitoring sites. 



Table 2. Values of NRIfocai, species richness, and abundances of Peucaea sparrows for each monitoring site. 



Peucaea sumichrasti Peucaea ruficauda 
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Site Richness NRIfocai NRIfocai Abundance NRI^cai NRIfocai Abundance 
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that the dissimilarity among communities was low when 
only taking species composition into account, but become 
highly dissimilar when abundance was taken into account. 



Thus, our findings are based on species turnover, which 
is only evident when it is analyzed through abundance. 
Individuals of one species interacted with individuals of 
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several species, and the level of competition differed in 
each interaction. We should think that individuals from 
one species are interacting with individuals of other spe- 
cies, rather than assuming that individuals interact with 
species; thus, abundances should be included in the analy- 
sis of phylogenetic structure whenever such data are avail- 
able (Norden et al. 2012).Therefore, competition intensity 
among organisms must be weighed by phylogenetic 
distance. Pragmatically, when the scale and data are 
appropriate, patterns and processes are, respectively, 
revealed and inferred by weighing phylogenetic distance 
by the number of individuals. 

The analyses were conducted using Peucaea population 
size as the species response; nonetheless, competition 
strength represented by the entire community could affect 
other population parameters, including birth, recruitment, 
and fecundity rates, among others. Gaston (2009) 
reviewed the population structure attributes necessary as 
a framework to continue research on the effects of com- 
munity strength on population range dynamics, and 
methods similar to the present study show promise. The 
phylogenetic community ecology approach as a means of 
studying a single species clarifies patterns and processes 
that otherwise might be masked. Although NRIfocai and 
NRI are related methodologies (r = 0.626, P < 0.05), the 
metrics told different stories about the same assemblages. 
The NRIfocai we applied here provided a new method to 
evaluate the potential community effect over the focal 
species. Both metrics are superficially similar and derived 
from the same theoretical framework; nonetheless, we 
caution that metric application and results differ and 
should be applied under different contexts. 
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Appendix Al: Accesion numbers 
included in the analysis. 


from GenBank for the species 


Peucaea ruficauda 


JQ1 73934 


Peucaea sumichrasti 


DQ433291 


Agelaius phoeniceus 


HM033218 


Amazilia Candida 


FJ027764 


Amazilia rutila 


EU442323 


Aratinga holochlora/strenua 


GU826182 


Burhinus bistriatus 


JQ1 74203 


Cacicus melanicterus 


JQ174231 


Calocitta formosa 


DQ433558 


Camptostoma imberbe 


DQ433420 


Campylorhynchus rufinucha 


DQ433425 


Carduelis psaltha 


JN801283 


Coccyzus minor 


JQ 174483 


Coiumbina inca 


DQ433B29 


Columba iivia 


JF498761 


Coiumbina passerina 


JN850709 


Coiumbina talpacoti 


JQ1 74508 


Coiinus virginianus 


DQ433524 


Corvus corax 


GUB71837 


Crotophaga suldrostris 


JN801306 


Cyanocompsa parellina 


FJ027473 


Cynanthus latirostris 


JN802021 
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Dives dives 


DQ433634 


Passerina caeruiea 


HIV1033638 


Diyocopus tinedtus 


JQ17472B 


Passer domesticus 


JQ1 75684 


Euplionia affinis 


EU442311 


Passerina iectiiancherii 


DQ433884 


Euphonia hirundinacea 


EU442315 


Patagioenas fiavirostris 


JN801347 


Geococcyx velox 


JN801317 


Piaya cayana 


JN801921 


icterus guiaris 


DQ433697 


Pitangus suipiiuratus 


JN801364 


icterus pectoraiis 


EU442293 


Poiioptiia aibiioris 


JN801366 


icterus pustuiatus 


DQ433688 


Pyrocepiiaius rubinus 


JQ288212 


Leptotila verreauxi 


FJ027748 


Quiscaius mexicanus 


EU442320 


Megarhynchus pitangua 


JN801790 


Sporophiiia minuta 


JQ1 76257 


Meianerpes aurifrons 


EU442305 


Sturneiia magna 


DQ433226 


Mimus giivus 


JN801810 


Tiiryotiiorus macuiipectus 


HIV1033838 


Mimus poiygiottos 


JQ1 75389 


Tiiryotiiorus pieurostictus 


HIV1208688 


Moiotiirus aeneus 


DQ433807 


Tityra semifasciata 


EU442296 


Moiottirus ater 


HIV1033B87 


Trogon citreoius 


JN802065 


Momotus mexicanus 


AY275849 


Turdus grayi 


JQ1 76566 


l\/lorococcyx eryttiropygus 


AY274054 


Turdus rufopaiiiatus 


HM033871 


Myiarchus nuttingi 


DQ433828 


Tyrannus meiancfioiicus 


JN802081 


l\/lyiozetetes simiiis 


FJ027897 


Voiatinia jacarina 


JQ627357 


Nyctidromus aibicoiiis 


JN801344 


Zenaida asiatica 


JQ176675 


Ortaiis poiiocepfiaia 


API 65496 







Appendix A2: Dissimilarity among monitoring sites using Bray-Curtis index. Upper diagonal reflects dissimilarity through qualitative data; 
meanwhile, lower diagonal is showing differences using abundances. 
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